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Abstract 

We propose a splitting method for solving equilibrium problems involving the sum of two 
bifunctions satisfying standard conditions. We prove that this problem is equivalent to find a 
zero of the sum of two appropriate maximally monotone operators under a suitable qualifica- 
tion condition. Our algorithm is a consequence of the Douglas-Rachford splitting applied to 
this auxiliary monotone inclusion. Connections between monotone inclusions and equilibrium 
problems are studied. 
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1. Introduction 

In the past years, several works have been devoted to the equilibrium problem 

find xeC such that (Vy G C) H(x,y)>0, (1.1) 

where C is a nonempty closed convex subset of the real Hilbert space H, and H : C x C — > H 
satisfies the following assumption. 

Assumption 1.1 The bifunction H: C x C — > R satisfies 
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(i) (Vz G C) = 0. 

(ii) (V(x, y) G C x C) y) + #(y, x) < 0. 

(iii) For every x in C, H(x, ■) : C — > R is lower semicontinuous and convex. 

(iv) (V(i,!/,z)6(7 3 ) rSn if((l -e)ar + £2,y) < H(x,y). 

e-»0+ 

Throughout this paper, the solution set of (1.1) will be denoted by Sh- 

Problem (1.1) models a wide variety of problems including complementarity problems, 
optimization problems, feasibility problems, Nash equilibrium problems, variational inequal- 
ities, and fixed point problems [1, 4, 5, 11, 16, 17, 18, 20, 24, 29]. Sometimes the bifunction 
H is difficult to manipulate but it can be considered as the sum of two simpler bifunctions 
F and G satisfying Assumption 1.1 (see, for example, [27]). This is the context in which we 
aim to solve problem (1.1). Our problem is formulated as follows. 

Problem 1.2 Let C be a nonempty closed convex subset of the real Hilbert space H. Sup- 
pose that F:CxC-fR and G:CxC-^l are two bifunctions satisfying Assumption 1.1. 
The problem is to 

find x G C such that (\/y G C) F(x,y) + G(x,y) > 0, (1.2) 
under the assumption that such a solution exists, or equivalently, Sf+g ^ - 

In the particular instance when G = 0, Problem 1.2 becomes (1.1) with H = F, which 
can be solved by the methods proposed in [16, 17, 20, 25, 28]. These methods are mostly 
inspired from the proximal fixed point algorithm [23, 34]. The method proposed in [33] can 
be applied to this case when F: (x,y) h- > (Bx \ y — x) and B is maximally monotone. On 
the other hand, when G: (x,y) h- > (Bx \ y — x), where B : H — > H is a cocoercive operator, 
weakly convergent splitting methods for solving Problem 1.2 are proposed in [11, 24]. Several 
methods for solving Problem 1.2 in particular instances of the bifunction G can be found 
in [6, 8, 9, 21, 30, 31, 32, 37] and the references therein. In the general case, sequential 
and parallel splitting methods are proposed in [26] with guaranteed ergodic convergence. 
A disadvantage of these methods is the involvement of vanishing parameters that leads to 
numerical instabilities, which make them of limited use in applications. The purpose of this 
paper is to address the general case by providing a non-ergodic weakly convergent algorithm 
which solves Problem 1.2. The proposed method is a consequence of the Douglas- Rachford 
splitting method [22, 36] applied to an auxiliary monotone inclusion involving an appropriate 
choice of maximally monotone operators. This choice of monotone operators allows us to 
deduce interesting relations between monotone equilibrium problems and monotone inclusions 
in Hilbert spaces. Some of these relations are deduced from related results in Banach spaces 
[2, 35]. 

The paper is organized as follows. In Section 2, we define an auxiliary monotone inclusion 
which is equivalent to Problem 1.2 under a suitable qualification condition, and some relations 
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between monotone inclusions and equilibrium problems are examined. In Section 3, we 
propose a variant of the Douglas-Rachford splitting studied in [3, 36] and we derive our 
method whose iterates converge weakly to a solution of Problem 1.2. We start with some 
notation and useful properties. 

Notation and preliminaries 

Throughout this paper, 7-L denotes a real Hilbert space, (■ | •) denotes its inner product, and 
|| • || denotes its induced norm. Let A: H — > 2 H be a set- valued operator. Then dom„4 = 
{x G % | Ax 7^ 0} is the domain of A and gra„4 = {(x, n) G H x H | u G Ax} is its graph. 
The operator A is monotone if 

(V(x, n) G gr& A)(\/(y, v ) G gra^4.) (x — y \ u — v) > 0, (1.3) 

and it is called maximally monotone if its graph is not properly contained in the graph of 
any other monotone operator in %. In this case, the resolvent of A, Ja = (Id +A)' 1 , is 
well defined, single valued, and dom^ = "H. The reflection operator IZj^ = IJjy — Id is 
nonexpansive. 

For a single-valued operator T: domT C M. — > H, the set of fixed points is 

Fix T = {x G H I x = Tx}. (1.4) 
We say that T is nonexpansive if 

(Vx G domT) (Vy G domT) \\Tx - Ty\\ < \\x - y\\ (1.5) 
and that T is firmly nonexpansive if 

(Vx G domT)(Vn G domT) ||Tx - Ty\\ 2 < \\x - y\\ 2 - ||(Id -T)x - (Id -T)y\\ 2 . (1.6) 



Lemma 1.3 (cf. [10, Lemma 5.1]) Lei T: domT = % — > H be a nonexpansive operator 
such that Fix 0. Let (/i„)„ e N a sequence in ]0, 1[ and (c n ) ne pj 6e a sequence in 7-L such 
that J2 neN fJ>n(l — Mn) = +00 and J2 n &N ^W^W < +°°- Let x G "H and set 

(Vn G N) x„ + i = x n + n n (Tx n + c n - x n ). (1.7) 
Then (x„)„ 6 pj converges weakly to x G Fix T and (x n — Tx„) ne N converges strongly to 0. 

Now let F: C x C — >■ R be a bifunction satisfying Assumption 1.1. The resolvent of T is the 
operator 

J F : H -> 2 C : x h-> {2 G C | (Vn G C) F(z,y) + (z - x | y - z) > 0} , (1.8) 
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which is single valued and firmly nonexpansive [11, Lemma 2.12], and the reflection operator 



R F : H ->■ H: x ^ 2.1 F x -x (1.9) 



is nonexpansive. 



Let C C H be nonempty, closed, and convex. We say that lies in the strong relative 
interior of C, in symbol, G sriC, if IJa>o ^ = spanC. The normal cone of C is the 
maximally monotone operator 

Afc:*^*:^!^!^' <»-*l«>^»}. if "^ (1 . 10 ) 

0, otherwise. 



We denote by r'o('H) the family of lower semicontinuous convex functions / from ~H to 
] — oo, +oo] which are proper in the sense that dom/ = {x e H \ f(x) < +00} is nonempty. 
The sub differential of / e To(H) is the maximally monotone operator df:'H^-2 H :xv-+ 
{u E 7i I (Vy 6K) (y — a; | u) + f(x) < f(y)}. For background on convex analysis, mono- 
tone operator theory, and equilibrium problems, the reader is referred to [3, 5, 11]. 



2. Monotone inclusions and equilibrium problems 

The basis of the method proposed in this paper for solving Problem 1.2 is that it can be 
formulated as finding a zero of the sum of two appropriate maximally monotone operators. 
In this section, we define this auxiliary monotone inclusion and, additionally, we study a class 
of monotone inclusions which can be formulated as an equilibrium problem. 

2.1. Monotone inclusion associated to equilibrium problems 

We first recall the maximal monotone operator associated to problem (1.1) and some 
related properties. The following result can be deduced from [2, Theorem 3.5] and [35, 
Proposition 4.2], which have been proved in Banach spaces. 

Proposition 2.1 Let F : C x C — > M be such that Assumption 1.1 holds and set 

A F :U^2 H :x^\{ UeU \ {yyeC) F &y) + (x-y\u)>0}, if x £ C; 

] 0, otherwise. 

Then the following hold: 

(i) Af is maximally monotone. 

(ii) Sp = zer Af- 
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(iii) For every 7 G ]0, +00 [, J^a f = J-yF- 



The following proposition allows us to formulate Problem 1.2 as an auxiliary monotone 
inclusion involving two maximally monotone operators obtained from Proposition 2.1. 

Theorem 2.2 Let C , F , and G be as in Problem 1.2. Then the following hold. 

(i) zei{A F + A G ) C S F+G . 

(ii) Suppose that span(C — C) is closed. Then, zer (Af + Ag) = Sf+g- 

Proof, (i). Let x G zer (Af + Ag)- Thus, x G C and there exists u G A F x D —Aqx, which 
yield, by (2.1), 

\WyeC) F(x,y) + (x-y\u)>0 
(\/y G C) G(x, y) + (x-y\-u)> 0. 
Hence, by adding both inequalities we obtain 

(yyeC) F(x,y) + G(x,y)>0 

and, therefore, x G Sp+ G - 

(ii). Let x G Sf+ g and define 



(2.2) 
(2.3) 



f:H—} ]— 00, +00] : y !->■ 



g: H — > ]— 00, +00] : y i-> 



F(x,y), 
-00, 

G(x,y), 
-00, 



if y e C; 
otherwise; 

if y e C; 

otherwise. 



(2.4) 



Assumption 1.1 asserts that / and g are in To('H), dom/ = dom(7 = C / 0, and since 
x G Sf+g, (1-2) yields / + g > 0. Hence, it follows from Assumption l.l(i) and (2.4) that 

0. (2.5) 



mm (f(y) + g{y)) = f(x) + g(x) 
yen v 



Thus, Fermat's rule [3, Theorem 16.2] yields G d(f + g)(x). Since span(C — C) is closed, 
we have G sri(C — C) — sri(dom/ — doing). Therefore, it follows from [3, Corollary 16.38] 
that G df(x) + dg(x) which implies that there exists u G H such that u G df(x) and 
—Mq G dg(x). This is equivalent to 



(VyeH) f(x) + (y-x\u )<f{y) 
(Wy G H) g(x) + (y - x | -w ) < g-(y). 

Since Assumption l.l(i) and (2.4) yield f(x) = g(x) = 0, we have that (2.6) is equivalent to 

'(WyeC) F(x,y) + (x-y\u o )>0 
(WyeC) G(x, y) + (x-y\ —u ) > 0. 

Hence, we conclude from (2.1) that u G Afx fl — Ag%, which yields x G zer(^l F + Ag)- U 



(2.6) 



(2.7) 
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2.2. Equilibrium problems associated to monotone inclusions 



We formulate some monotone inclusions as equilibrium problems by defining a bifunction 
associated to a class of maximally monotone operators. In the following proposition we 
present this bifunction and its properties. 

Proposition 2.3 (c/. [11, Lemma 2.15]) Let A: H — > 2 n be a maximally monotone operator 
and suppose that C C int dom„4. Set 



Then the following hold: 

(i) Fa satisfy Assumption 1.1. 

(ii) J Fa = Ja+Mc- 

Remark 2.4 Note that the condition C C int dom .4. allows us to take the maximum in 
(2.8) instead of the supremum. This is a consequence of the weakly compactness of the sets 
(Ax) xe c ( see [11, Lemma 2.15] for details). 

Proposition 2.5 Let A: H — > 2 H be a maximally monotone operator and suppose that C C 
int dom A. Then zer (A + Mc) = Sf a ■ 

Proof. Indeed, it follows from [3, Proposition 23.38], Proposition 2.3, and [11, Lemma 2.15(i)] 
that 



which yields the result. □ 
Remark 2.6 

(i) Note that, in the particular case when dom A = int dom A = C = H, Proposition 2.5 
asserts that zer A = Sp A , which is a well known result (e.g., see [19, Section 2.1.3]). 

(ii) In Banach spaces, the case when C = dom A C H is studied in [2, Theorem 3.8]. 

The following propositions provide a relation between the operators defined in Proposi- 
tions 2.1 and 2.3. 



Proposition 2.7 Let B: % — > 2 H be maximally monotone and suppose that C C int dom B. 
Then, A Fb = B + M c . 




: (x, y) i — y max (y — x \ u). 

u£Ax 



(2.8) 



zer (.4 + M c ) = Fix (J A+A r c ) = Fix ( J Fa ) = S Fa , 



(2.9) 
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Proof. Let (x,u) G "H 2 . It follows from (2.8) and [5, Lemma 1] (see also [11, Lemma 2.14]) 
that 



u G Af b x x G C and (Vy G C) F B (x, y) + (x — y | w) > 

<^> x G C and (\/y G C) max (y — x \ v) + (x — y | u) > (2.10) 

vGBx 

xeC and (Vy G C) max (y - a; | v - u) > ( 2 -H) 

^ iGC and (3v G i3x)(Vy G C) (y - x | f - u) > (2.12) 
ueBx+Ncx, (2.13) 

which yields the result. □ 

Proposition 2.8 LetG be such that Assumption 1.1 holds, and suppose that C = dom^c = 
H. Then, F Ag <G. 

Proof. Let (x,y) G CxC and let u G Aqx. It follows from (2.1) that G(x, y) + (x — y \ u) > 0, 
which yields 

(Vm G (j/ - ar I u) < G(x, y). (2.14) 

Since C = int dom^4 G = H, the result follows by taking the maximum in the left side of the 
inequality. □ 

Remark 2.9 

(i) Note that the equality in Proposition 2.8 does not hold in general. Indeed, let % = M, 
C = H, and G: (x,y) (-> y 2 — x 2 . It follows from [11, Lemma 2.15(v)] that G satisfy 
Assumption 1.1. We have u G Aqx (Wy G H) y 2 — x 2 + (x — y \ u) > -v^ u = 2x 
and, hence, for every (x, y) G % x %, F^ G (x, y) — (y — x)2x = 2xy — 2x 2 . In particular, 
for every y G R\ {0}, F Ag (0, V ) = < y 2 = G(0,y). 

(ii) In the general case when C = dom^4 C H, necessary and sufficient conditions for the 
equality in Proposition 2.8 are provided in [2, Theorem 4.5]. 



3. Algorithm and convergence 

Theorem 2.2(h) characterizes the solutions to Problem 1.2 as the zeros of the sum of two 
maximally monotone operators. Our algorithm is derived from the Douglas-Rachford split- 
ting method for solving this auxiliary monotone inclusion. This algorithm was first proposed 
in [13] in finite dimensional spaces when the operators are linear and the generalization to 
general maximally monotone operators in Hilbert spaces was first developed in [22]. Other 
versions involving computational errors of the resolvents can be found in [10, 14]. The conver- 
gence of these methods needs the maximal monotonicity of the sum of the operators involved, 
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which is not evident to verify [3, Section 24.1]. Furthermore, the iterates in these cases do 
not converge to a solution but to a point from which we can calculate a solution. These 
problems were overcame in [36] and, later, in [3, Theorem 25.6], where the convergence of the 
sequences generated by the proposed methods to a zero of the sum of two set- valued operators 
is guaranteed by only assuming the maximal monotonicity of each operator. However, in [36] 
the errors considered do not come from inaccuracies on the computation of the resolvent but 
only from imprecisions in a monotone inclusion, which sometimes could be not manipulable. 
On the other hand, in [3, Theorem 25.6] the method includes an additional relaxation step 
but it does not consider inaccuracies in its implementation. 

We present a variant of the methods presented in [36] and [3, Theorem 25.6], which 
has interest in its own right. The same convergence results are obtained by considering a 
relaxation step as in [14, 15] and errors in the computation of the resolvents as in [10, 14]. 

Theorem 3.1 Let A and B be two maximally monotone operators from rt to 2 n such that 
zer(.4 + B) 7^ 0. Let 7 G ]0,+oo[ ; let (\ n ) ne n be a sequence in ]0,2[, and let (a n ) n6 N and 
(frn)neN be sequences in H such that b n — 0, 

2J A n (2 - A n ) = +00, and 2J A n (||a n || + ||6 n ||) < +00. (3.1) 

neN ngN 

Let Xq Erl and set 

Zn = JjA^Vn ~ x n ) + a n (3.2) 
Xn+1 X n "I - Aj 1 (^ ri Vn) ■ 

Then there exists x G Fix {TZ 1 j{IZ 1 b) such that the following hold: 

(i) J lB x G zer(„4 + B). 

(ii) {Ti^j^ilZ^BXn) — i n )neN converges strongly to 0. 
(hi) (x„) ne N converges weakly to x. 

(iv) (y n ) ne N converges weakly to J^b^- 

Proof. Denote T = 'R-^A^-yB- Since TZ^a and IZ^b are nonexpansive operators, T is nonex- 
pansive as well. Moreover, since [3, Proposition 25.1(h)] states that J^Fix T) = zer(.A+£>), 
we deduce that Fix 0. Note that (3.2) can be rewritten as 

(Vn G N) x n+ i = x n + fi n (Tx n + c n - x n ), (3.3) 

where, for every n G N, 

= ^ and c n = 2(j jA (2(J^ B Xn + b n ) - x n ) + a n - J lA (2J y BX n - x n ) ~ b n ). (3.4) 
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Hence, it follows from the nonexpansivity of J lA that, for every n G N, 

\\Cn\\ = 2\\J lA (2(J lB X n + b n ) - X n ) +d n - J^A^JjBXn - x n ) - b n \\ 

< 2\\J lA (2(J~ fB x n + b n ) - x n ) - J-yA^J-rBXn - x n )\\ +2\\a n \\+ 2\\b n \\ 

< 2\\2(J lB x n + b n ) -x n - (2J jB x n - x n )\\ + 2||a n || + 2||6„|| 

= 2{\\a n \\ +3||6 n ||) (3.5) 

and, therefore, from (3.1) and (3.4) we obtain 

2jjMe n || < z2^n(\\an\\ +3||6 n ||) < 3^jA n (||a n || + ||6 n ||) < +oo. (3.6) 

neN nGN neN 

Moreover, since the sequence (A n ) ng N is in ]0,2[, it follows from (3.4) that (/i n )neN is a 
sequence in ]0, 1[ and, from (3.1) we obtain 

J2vn(l-Vn) = ^A n (2-A n ) = +oo. (3.7) 

neN neN 

(i) . This follows from [3, Proposition 25.1(h)]. 

(ii) and (iii). These follow from Lemma 1.3. 
(iv). From the nonexpansivity of J^ B we obtain 

\\Vn - 2/0 1| < \\JjBX n ~ JjBXo\\ + ||&n ~ M < \\x n ~ X \\ + ||6„ - 6 ||. (3.8) 

It follows from (iii) and b n — ^ that (x n ) ne ^ and (6 n ) ne N ar e bounded, respectively. Hence, 
(3.8) implies that (y n ) n eN is bounded as well. Let y G H be a weak sequential cluster point 

of (y n )nen, say y kn -± y, and set 



(Vn G N) < 



£/n J^jBXn 

%n yJ^jS^Un Xn) 

l^n 2?/ n X n Z n 

x n y n . 



(3.9) 



It follows from (3.2) that 



{z n ,u n ) G gra7^l 
(V//GN) ^ (y n ,v n ) G gra7^ (3.10) 

"I - ^n 2/n ^n- 
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For every n G N, we obtain from (3.9) 

\\Zk n -yk n \\ = \\JjA(2J-yBXk„ — Xkn) - JjBXknW 

— 2 W^jA^^B^kn — %kn) ~ ( 2 J~ 7 B£fc n — Xk„) - Xk n \\ 

= ^P^A{^B^k n ) -Xkn\\. (3.11) 

Hence, (ii) yields z kn — Vk n — > 0, and, therefore, from (3.10) we obtain that Uk n + v kn — > 0. 
Moreover, it follows from 6 fcn — 0, y\.„ V, an d (3.2) that y kn — y, and, hence, Zk n — ^ y. 
Thus, from (iii) and (3.9), we obtain Uk n — y — x and Vk n — x — y. Altogether, from [3, 
Corollary 25.5] we deduce that y G zer(7^4 + jB) = zer(.A + B), (y,y — x) G gra7-4, and 
(y, x — y) G gra7$. Hence, y = J^bx and y G domA Therefore, we conclude that J^bx is 
the unique weak sequential cluster point of (y n )nm and then y n — ^ J^bx. □ 

Now we present our method for solving Problem 1.2, which is an application of Theorem 3.1 
to the auxiliary monotone inclusion obtained in Theorem 2.2. 

Theorem 3.2 Let C , F , and G be as in Problem 1.2 and suppose that span(C — C) is closed. 
Let 7 G ]0,+oo[, let (\ n ) n &i be a sequence in ]0,2[, and let (a„) ne ^ and (6 n ) ne N be sequences 
in 7-L such that b n 0, 

^A n (2-A„) = +oo, and ^A„(||a n || + ||6 n ||) < +oo. (3.12) 

neN neN 
Let Xq G H and set 

{Vn = J-yGXn + b n 
z n = J lF (2y n - x n ) + a n (3.13) 
Xn+l X n -\- \ n (^Z n yn) ■ 

Then there exists x G Fix (RjfRjg) such that the following hold: 

(i) J^qx G Sf+g- 

(ii) (RypffijoXn) ~ x n ) n< zfq converges strongly to 0. 

(iii) (x n ) n£ N converges weakly to x. 

(iv) (y n ) n m converges weakly to J 7 gx. 

Proof. Note that, from Theorem 2.2(h), we have that 

zer(^l F + Ag) = S F+G ^ 0, (3. 14) 

where Af and Ac are defined in (2.1) and maximally monotone by Proposition 2.1 (i) . In 
addition, it follows from Proposition 2.1 (iii) that (3.13) can be written equivalently as (3.2) 
with A = Af and B = Ag- Hence, the results are derived from Theorem 3.1, Proposition 2.1, 
and Theorem 2.2. □ 
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Remark 3.3 Note that the closeness of span(C — C) and Theorem 2.2(h) yields (3.14), which 
allows us to apply Theorem 3.1 for obtaining our result. However, it is well known that this 
qualification condition does not always hold in infinite dimensional spaces. In such cases, it 
follows from Theorem 2.2(i) that Theorem 3.2 still holds if zer(AF + Ag) ^ 0- Conditions 
for assuring existence of solutions to monotone inclusions can be found in [7, Proposition 3.2] 
and [3]. 

Finally, let us show an application of Theorem 3.2 for solving mixed equilibrium problems. 
Let / G r ("H). For every x G 71, prox^x is the unique minimizer of the strongly convex 
function y \-¥ f(y) + \\y — x\\ 2 /2. The operator prox^: % — > H thus defined is called the 
proximity operator. 

Example 3.4 In Problem 1.2, suppose that G: (x,y) h-> f(y) — f(x), where / G To^H) is 
such that C C dom/. Then Problem 1.2 becomes 

find x G C such that (Vy G C) F(x,y) + f(y) > f(x), (3.15) 

which is known as a mixed equilibrium problem. This problem arises in several applied 
problems and it can be solved by using some methods developed in [8, 32, 30, 37]. However, 
all this methods consider implicit steps involving simultaneously F and /, which is not easy 
to compute in general. On the other hand, it follows from [11, Lemma 2.15(v)] that (3.13) 
becomes 

{y n = prox tc+7/ x n + b n 
z n = J lF {2y n - x n ) + a n (3.16) 
Xn+l X n ~\- X n ( K Z n 2/n); 

which computes separately the resolvent of F and the proximity operator of /. If span(C — C) 
is closed, Theorem 3.2 ensures the weak convergence of the iterates of this method to a 
solution to (3.15). Examples of computable proximity operators and resolvents of bifunctions 
can be found in [12] and [11], respectively. 
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